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Abstract 



We explore the quark-mass dependence of two-nucleon systems. Allowed 
regions for the scattering lengths in the and 3 S± channels as functions of 
the light-quark masses are determined from the current uncertainty in strong- 
interaction parameters that appear at next-to-leading order in the effective 
field theory. Where experimental constraints are absent, as is the case for 
the quark-mass dependent four-nucleon operators, we use naive dimensional 
analysis. We find it likely that there is no bound state in the channel in 
the chiral limit. However, given the present uncertainties in strong-interaction 
parameters it is unclear whether the deuteron is bound or unbound in the 
chiral limit. 



I. INTRODUCTION 



An important aspect of strong-interaction physics that is only now being explored is the 
dependence of nuclei and nuclear reactions on the fundamental parameters of nature, such 
as Aq CD and the quark masses, m q . With recent developments in the effective field theory 
(EFT) description of multi-nucleon systems fundamental questions about nuclear physics 
can be addressed. While the m g -dependence of the nuclear force is unrelated to present 
day observables, it is a fundamental aspect of nuclear physics, and in some sense serves as 
a benchmark for the development of a perturbative theory of nuclear forces. Having this 
behavior under control will be essential to any bridge between lattice QCD simulations and 
nuclear physics in the near future. In work by Bedaque, van Kolck and the authors [1] 
(which we will refer to as BBSvK), a new power-counting scheme, BBSvK-counting, was 
introduced, and with a partial next-to-leading order (NLO) calculation in this new counting 
scheme, an attempt was made to determine the behavior of the two-nucleon force as the m q 
were increased or decreased from their physical values. Calculations of the m g -dependence 
of scattering in the ^-channel and 3 Si — 3 D 1 coupled-channels as well as a discussion of the 
deuteron binding energy were presented. In a subsequent paper [2] we performed the first 
complete NLO calculation in BBSvK-counting of the m^-dependence of the two-nucleon sys- 
tems. Previous estimates of such m g -dependences [3,4] that served as input to place bounds 
on physics beyond the standard model via the time-variation of fundamental parameters on 
cosmological time scales were found to be incomplete. 

In this work, which should be considered a sequel to Ref. [2], we investigate the m q - 
dependence of the two-nucleon systems, incorporating the present uncertainties in strong 
interaction parameters *. In some instances, parameters that play a central role in the 
m^-dependence, such as the leading m g -dependent four-nucleon interaction, are not con- 
strained experimentally. In such cases, we use naive-dimensional-analysis (NDA) to provide 
a "reasonable range" for their values. Of particular interest to us is the behavior of the 
two-nucleon sector in the chiral limit. We find it likely that there is no di-nucleon bound 
state in the ^o-channel; i.e. the scattering length remains negative. On the other hand, we 
cannot determine if there is a bound-state in the 3 S 1 — 3 £>i coupled channels in the chiral 
limit. The strong interaction parameters are sufficiently uncertain at present to preclude a 
definitive answer to this fundamental question. 

II. THE %-CHANNEL AND THE DI-NUCLEON 

In the ^-channel, pions give a subleading contribution to the scattering amplitude 
and are included in perturbation theory. In this channel, BBSvK-counting is equivalent to 
Kaplan-Savage- Wise (KSW) power-counting [5]. Thus, analytic expressions for the scatter- 
ing amplitude, and hence the scattering length, as a function of m q and momentum, p, are 
straightforwardly found. From the NLO amplitude [5] it is easy to construct pcot8^ s °\ 



*We are grateful to J. Bjorken who suggested this work to us during Marshall Baker Symposium 
held at the UW in June 2002. 
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which has a well-behaved power-series expansion for p < m n /2, and thus a linear com- 
bination of Co (the coefficient of the m q - and p-independent four-nucleon operator) and 
D 2 (the coefficient of the p-independent and leading m g -dependent four-nucleon operators) 
can be determined in terms of the scattering length a^ ^ = —23.714 ± 0.013 fm at the 
physical value of the pion mass. Furthermore, C 2 (the coefficient of the m^-independent, 
leading p 2 -dependent four-nucleon operator) can be determined in terms of the effective 
range, = 2.73 ± 0.03 fm. Once these parameters are fixed, the scattering length, ef- 

fective range and phase-shift can be determined as a function of m q at NLO in the KSW 
expansion. The scattering length is, in the isospin limit, 
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where 7 = p + 4n / 'M^Co(p), and p is the PDS renormalization scale [5]. We have used 
the leading order (LO) relation between the pion mass, m n , and the light-quark masses, 
m q . As pions are treated in perturbation theory, the m^-dependence arising from the chiral 
expansions of qa, f n an d M N is formally higher order in KSW-counting and is neglected. 
The physical values of these parameters, qa ~ 1-25, f w ~ 132 MeV and M N ~ 939 MeV, 
are used for the analysis in this channel. 

At present, the coefficient D 2 has not been experimentally separated from the coefficient 
Cq in this channel, or in the 3 S± — 3 Di coupled-channels. This is because both parameters 
are coefficients of momentum-independent operators with the only difference between the 
operators arising from their couplings to pions; the Co operator does not give rise to inter- 
actions between two-nucleons and multiple pions, while the D 2 operator does. In order to 
determine the range of values of D 2 consistent with NDA, we use the expression for the scat- 
tering length in eq. (1) at the physical value of the pion mass and find values of C (m™ YS ) 
and _D2( m ^ HYS ) that satisfy the constraint 
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where the constant r\ is chosen to be 1/5 and 1/15 for demonstrative purposes. It is worth- 
while discussing why these values are chosen. Very roughly, r\ ~ (m w /A) 2 where A represents 
the scale at which the EFT description breaks down. Naively, the KSW perturbative-pion 
expansion breaks down at a scale Ann = ^fn/dA^ ~ 300 MeV. This would then give 
i] ~ 1/5. On the other hand, one might expect that a better estimate of the cutoff scale 
is given by the mass of the fictitious scalar, a, that provides intermediate-range attraction 
in many potential models. Taking m a ~ 500 — 600 MeV one then arrives at rj ~ 1/15. Of 
course this last argument is somewhat specious; it is really at m a /2 ~ Ann that a should 
be treated as a dynamical object. This discussion underscores the difficulty in providing an 
a priori estimate of a parameter in an EFT. Hence, given the complete lack of experimental 
information on D 2 , and the absence of any reliable method to estimate its value, we take r\ 
to be 1/5 and 1/15 to provide an estimate of the uncertainty arising from D 2 . We stress that 
even larger values of 77, such as 1/3, are not excluded experimentally. Of course larger values 
of rj lead to a smaller radius of convergence of the EFT. In Fig. 1 we show the scattering 
length in the ^-channel as a function of m g , where we have randomly sampled values of 
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FIG. 1. The scattering length in the ^o-channel as a function of the pion mass. The light gray 
region corresponds to r/ = 1/5 and the black region corresponds to rj = 1/15. At the physical value 
of the pion mass the scattering length is a^ So ^ = —23.714 ± 0.013 fm. 

and 77 within the allowed ranges^. From the standpoint of NDA, it is likely that the 
di-neutron is unbound in the chiral limit for any values of m q less than the physical values. 
Of course, the same will be true for the other di-nucleon states in the ^-channel by isospin 
symmetry. On the other hand, for m q larger than their physical values, we see that there 
are both bound and unbound di-neutrons for values of the strong-interaction parameters 
consistent with NDA. This observation reinforces our previous analysis [2]. 

For the strong-interaction parameter sets where the scattering length remains unnaturally 
large in the chiral limit, it is likely that higher-order contributions to the scattering amplitude 
will significantly disturb the system. Thus, for such situations, higher-order contributions 
may lead to a bound state where one does not exist at NLO. 

III. THE 3 Si - 3 £>i COUPLED-CHANNELS AND THE DEUTERON 

The deuteron resides in the 3 S\ — 3 Di coupled-channels when the m q assume their physical 
values. In our previous work [1,2] we found that insufficient is presently known about the 
strong-interaction couplings in the two-nucleon sector to determine if the deuteron is bound 
or unbound for small values of m q , including the chiral limit, or for m q larger than their 
physical values. In this section we make these findings more explicit and determine the 
regions allowed by data, and by NDA in the absence of data, for the 3 Si scattering length 
and deuteron binding energy, B^, as a function of m q , at NLO in BBSvK-counting. 

The details of the calculation of scattering lengths, phase shifts and bound state energies 
in the 3 Si — 3 Di coupled channels can be found in Refs. [1,2], and we do not repeat the details 



^We choose to use scatter plots as the point density reflects the probability associated with a 
particular set of low-energy constants. 



here. However, a brief description of the procedure, and of the sources of uncertainty is in 
order. In this channel, unlike the ^-channel, the chiral limit of one-pion-exchange (OPE) 
is included at LO in BBSvK-counting. This then allows for the re-summation of the large 
tensor-force that persists in the chiral limit, and is responsible for the non-convergence of 
KSW-counting in this channel [6]. Deviations from the chiral limit are included in pertur- 
bation theory, as are higher-dimension operators and multiple pion exchanges [1]. In the 
3 S± — 3 Di coupled-channels, OPE generates both central and tensor potentials, 

Vc ' ( r ? m ?r) = ~ a * m l » V^\r]m^) = -a n f — + + m 2 ) , (3) 

where a w = g\(l — 2m^d 18 / g A ) 2 / (8nf 2 ). The low-energy constant (LEC) d 18 has been 
determined through various fitting procedures [7-10] to be —0.78 ± 0.27, —0.83 ± 0.06, 
-1.30 ± 0.24 [8,9] and -10.14 ± 0.45 GeV~ 2 [10] at third-order in the chiral expansion. 
While we are unable to establish a definition of the uncertainties associated with each of 
these determinations, we assume that they represent la errors arising from a parabolic 
error analysis of the MINUIT package, as partially described in Ref. [10]. Further, in the 
absence of an error correlation matrix, we will treat the uncertainties in the values of LECs 
as uncorrelated; this constitutes a weakness in our propagation of errors. Thus, we sample 
values of rf 18 over the range —1.54 GeV~ 2 < d 18 < —0.51 GeV -2 , corresponding to the la 
limits of the first three determinations above *, and as discussed in Ref. [9], this range is 
probably an underestimate of the true uncertainty in d 18 . A fourth-order calculation has been 
performed and extracted values of di$ that are consistent with the three values above [9]. 

As the chiral limit of the potentials in eq. (3) contribute at LO, the leading m q - 
dependence of g A , fn and Mn are required at NLO, as discussed in BBSvK. Each of these 
observables has been studied extensively, the results of which can be found in Refs. [7,11-13], 
and up to next-to-next-to-leading order (NNLO) it is known that 
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where U = 4.4 ± 0.2 [11,12], c\ ~ —1 GeV 1 [7] are m^-independent constants (we have 
explicitly separated the logarithmic contribution from 7 4 ) and the superscript (0) denotes 
the chiral limit value. We use g A = 1.25, M N = (M n + M p )/2 and f w = 135 MeV. The 
parameter in eq. (4) that is most uncertain and gives rise to a substantial uncertainty 
in our calculations in the two-nucleon sector is die- A complete analysis by Fettes [9] of 
the ttN sector provides three different determinations of die'. —0.91 ± 0.74, —1.01 ± 0.72 
and —1.76 ± 0.85 GeV" 2 . Thus we sample this parameter over the range —2.61 GeV -2 < 



*The value, d 18 10 GeV" 2 [10] is significantly larger than the other determinations and we 

discard it. 
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die < —0.17, again corresponding to the la limits of the three determinations. In our 
previous paper [2], we retained only the chiral logarithm contribution to g^, and evaluated 
its argument at the scale A = 500 MeV, which corresponds to choosing die ~ +1 GeV~ 2 . 
Note that while this value is outside the la range of g?i 6 , it is not excluded by data. As I 4 
makes only a small contribution to the OPE potential, its uncertainty has negligible impact 
and we hold I4 fixed to its central value in what follows. 

At NLO in BBSvK-counting there is a contribution from the chiral limit of two-pion 
exchange (TPE) and from an insertion of the C 2 and D 2 operators. The TPE potential in 
coordinate space has been computed in Ref. [14-16], and in the chiral limit becomes 

vtV, ) = m9 \l l :;r l) ± , ^>(r;0) = -^i. (5) 

In order to regulate the singular potentials that occur at LO and NLO we use a spatial 
square-well of radius R [1,2,17], where the potential outside the square well is 

/. s_( -M N V c (r\m v ) -2^/2M N V T ( r ', m Tr) \ (n\ 

L{r] ~ { -2V2 M N V T (r; m n ) -M N (V c (r; m n ) - 2V T (r; m w )) - Q/r* J ' W 

where 

V c (r;m K ) = V§'\r;m K ) + V^ } (r;0) , V T (r;m w ) = V^(r;m n ) + ^(^0) . (7) 
The potential inside the square well is 

v /_. _ k 2,_( -M N (V Co + V D2 ml + k 2 V C2 ) \ 

where Vc , Vd 2 and Vc 2 are constant potentials corresponding to the renormalized local 
operators with coefficients Co, D 2 and C2 in the 3 <Si — 3 Di coupled-channels, respectively 
and again we have used the LO relation between the pion mass and m q . As discussed in 
BBSvK, we can make an identification between the coefficients of local operators and the 
constant potentials of the square- wells that enter into eq. (8), e.g. 

Q^) (r) ^ 3Q^-r) ^ Vct6{R _ r) (g) 

and similarly for Dj and Vo v It is important to recall that there is implicit m 9 -dependence 
in this potential due to qa, and f n , in addition to the explicit dependence from D 2 m^, 
and from OPE. Defining the wavefunction \& to be 



§The appearance of momentum-independent step-functions in the D-wave potential might appear 
somewhat puzzling. However, this arises naturally when a spatial regulator is used to "smear out" 
the local short-distance part of the central potential in a manner that is consistent with rotational 
invariance. Contributions of this D-wave component to the scattering amplitude are suppressed 
by powers of the square- well radius, R, as R — > 0, as required in the renormalized theory. 
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where u(r) is the S-wave wavefunction and w(r) is the D-wave wavefunction, the regulated 
Schrodinger equation is 

*"(r) + [k 2 + V L (r;m n ) 6(r - R) + V s (r;m w , k 2 ) 9(R - r) ]^(r) = . (11) 

It is this Schrodinger equation that we solve to generate the scattering length and location of 
the bound-state (if present) in this channel. A detailed comparison of phase shifts generated 
with this formulation, and the results of the Nijmegen partial wave analysis [18] can be 
found in Ref. [1,2]. 

A few words are in order concerning the power-counting. BBSvK-counting concerns itself 
with a consistent removal of cutoff effects, order- by-order in the EFT expansion [1]. The 
OPE potential is special in that only the chiral limit piece requires renormalization. This 
is because the 1/r 2 contribution in the pion mass expansion vanishes and all other terms in 
the expansion are regular at the origin. The chiral expansion of OPE has been shown to 
converge, albeit slowly, in Ref. [1]. We could, of course, arrange the power counting so as 
to keep as many orders in the OPE perturbative chiral expansion as necessary in order to 
match to the results shown in our paper, which treats the chiral expansion of OPE to all 
orders; however, this would be rather pointless, since treating the potential to all orders in 
the chiral expansion is perfectly innocuous from the perspective of renormalization, and in 
such cases it is perfectly legitimate to treat perturbative physics nonperturbatively. On the 
other hand, the TPE potential has a significantly more divergent chiral expansion as there 
are singular-attractive 1/r 5 and 1/r 3 contributions. It is likely that treating TPE to all 
orders while keeping only the counterterms at NLO in Weinberg counting does not remove 
all cutoff artifacts (scaling as logarithms or inverse powers of the coordinate-space cutoff). 
We have confirmed this strong cutoff dependence of the full TPE numerically; this is, in 
fact, strong motivation for using BBSvK counting. 

Using the same definition of rj as in eq. (2), we generate values for D 2 and Co consistent 
with NDA, and solve the Schrodinger equation in eq. (11) for the scattering length in the 3 S± 
channel and the deuteron binding energy. We choose a cutoff of R = 0.45 fm. The cutoff- 
independence of our results is discussed in Refs. [1,2]. In Fig. 2, the scattering length in the 
3 Si channel is shown as a function of the pion mass, where we have used the LO relation 
between m q and m n . Again we have chosen 77 to have the values r] — 1/5 and 1/15 to illustrate 
the effect of D 2 . We take d 18 = -1.54 GeV~ 2 and -2.61 GeV" 2 < d m < -0.17 GeV~ 2 [9]. 
We again randomly and uniformly sample parameter space within the allowed range**. In 
Fig. 3 the scattering length in the 3 S 1 channel is shown as a function of the pion mass. The 
difference with Fig. 2 is that in Fig. 3 we take d 18 = —0.51 GeV~ 2 In Fig. 4 we show 



**A more appropriate procedure might be to sample d±Q from a Gaussian-weighted distribution. 

ttWhile dis does not contribute in the chiral limit, it does alter the strength of the nucleon-nucleon 
force at the physical value of the pion mass, where the local four-nucleon operators are fit to data. 
Data, of course, contains all orders in the chiral expansion. 
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FIG. 2. The scattering length in the 3 Si-channel as a function of the pion mass. The light 
shaded region corresponds to 77 = 1/5 and the black region corresponds to rj = 1/15. The param- 
eter d\e, defined in eq. (4), is taken to be in the interval —2.61 GeV -2 < die < —0.17 GeV -2 [9], 
and dig = —1.54 GeV -2 . At the physical value of the pion mass the scattering length is 
a( 3s i) ~ +5.425 fm. 




FIG. 3. The scattering length in the 3 Si-channel as a function of the pion mass. Parameters 
and labels are the same as Fig. 2 except dig = —0.51 GeV -2 . 
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FIG. 4. The deuteron binding energy as a function of the pion mass (from the physical value to 
the chiral limit). We do not show values below = 60 MeV since the deuteron can be both bound 
and unbound. The shaded region corresponds to rj = 1/5, —2.61 GeV -2 < die < —0.17 GeV~ 2 , 
and -1.54 GeV~ 2 < dig < -0.51 GeV" 2 . 
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FIG. 5. The scattering length in the 3 Si-channel as a function of the pion mass for 
d w = +1 GeV~ 2 , r] = 1/3 and d i8 = -0.51 GeV" 2 . 
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the deuteron binding energy as a function of the pion mass for the full range of d\ 8 and d 16 
and taking 77 = 1/5. 

Perhaps it is worth interpreting the binding energy curves in Fig. 4 using the scattering 
length curves of Fig. 2 and Fig. 3. It is obvious from Fig. 2 and Fig. 3 that the scattering 
length can be negative for m n < 60 MeV, indicating an unbound deuteron. Naively one 
would expect the scattering length to go through infinity to unbind the deuteron; this would 
correspond to the deuteron binding energy going through zero. However, in this work we 
find that the deuteron becomes unbound in the approach to the chiral limit by becoming 
increasingly deeply bound and thus leaving the range of applicability of the EFT. 

From these figures we conclude that the present uncertainty in strong-interaction param- 
eters precludes a definitive statement being made concerning the binding or non-binding of 
the deuteron in the chiral limit. This is consistent with our previous studies [1,2]. It is 
a combination of uncertainties in the 7rN sector and the lack of any information in the 
two-nucleon sector regarding D2 that conspire to prevent a definitive conclusion. It is clear 
that progress must be made on both issues before the chiral limit of the 3 Si — Z D\ coupled 
channels can be predicted. To demonstrate the impact of die, in Fig. 5 we show the scat- 
tering length in the 3 Si channel for d 16 = +1 GeV -2 , a value that is not excluded by data. 
Clearly the sensitivity of the two-nucleon system to this parameter is rather dramatic. For 
this parameter set we see the deuteron unbinding by the scattering length going through 
both infinity and zero. 

Impressive progress has been made in extrapolating lattice observables in the single- 
nucleon and meson sectors from lattice quark masses to the physical quark masses. In 
addition to being important in its own right, this work is essential for interpretation of 
lattice QCD results in nuclei. Extrapolations of some quenched lattice QCD observables 
at large quark masses (m^ ~ 500 MeV) to physical values of the quark masses are seen 
to be in good agreement with experimental values [19]. Future partially-quenched or un- 
quenched simulations are required at smaller quark masses in order to have confidence in 
first principles calculations. As we discussed in Ref. [1,2], a quenched lattice calculation 
of the scattering lengths in both the ^So and ^i-channels exists [20], but for pion masses 
greater than ~ 550 MeV. While it would be tempting to extrapolate our results to these 
pion masses it is likely that the EFT will have broken down at much smaller pion masses. 
Moreover, as pointed out in Ref. [21], the long-distance part of the two-nucleon force is mod- 
ified in quenched and partially-quenched QCD; for instance, long-distance Yukawa behavior 
is modified to pure exponential fall-off. This poses yet another challenge to bridging the gap 
between nuclear physics and lattice QCD [21]. 

IV. COMMENTS ON RECENT WORK BY EPELBAUM, 
GLOCKLE AND MEIfiNER 

While this paper was in preparation, an analysis by Epelbaum, Glockle and MeiBner 
(EGM) has become available [22] which also extends earlier work on this subject [1,2]. 
EGM consider the m^-dependence of nuclear forces using their formulation of Weinberg's 
power-counting [23]. In contrast with our coordinate-space derivation, this work involves a 
numerical solution to the Lippmann-Schwinger equation with a potential generated order-by- 
order in the chiral expansion with a power-counting derived for the meson and single nucleon 
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sectors. Their work is in qualitative agreement with the results of our previous papers [1,2] 
and with the results we have presented here. Despite the differences in regularization method 
and the subtly different power-counting schemes, it not surprising to see that the two EFTs 
yield results that appear to be perturbatively close. In the ^-channel, it is encouraging 
to see that the results of the numerical calculations of Ref. [22] agree with our analytic 
results obtained from the perturbative-pion expansion. However, there are quite significant 
differences in the uncertainties associated with the extrapolation of the observables away 
from the physical value of the pion mass, particularly in the 3 <Si — 3 Di coupled channels. 
These discrepancies arise from a different treatment of uncertainties associated with the 
strong-interaction parameters d 16: d 18 and D 2 . It appears that EGM [22] have chosen a 
range of values of Di consistent with r\% 1/19, which we consider to be overly optimistic 
EGM argue that, for the range of cutoffs they use to define their EFT, all experimentally 
determined low-energy constants in the two-nucleon sector fall within the "NDA" estimates 
that they have created. Unfortunately, at present there is no way to test the reliability of 
such estimates for m^-dependent operators. Moreover, EGM have chosen not to propagate 
the known and relatively large errors associated with the extraction of d\§ [9] from the single 
nucleon sector. In Fig. 6 we show the 3 S± scattering length as a function of the pion mass 
for the parameter ranges used by EGM. With their parameters we reproduce their results, 
up to terms higher order in both EFT descriptions. For instance, we recover their result of 
1.25 fm< a( 35l )< 1.85 fm in the chiral limit. Therefore, while we find the claims of Ref. [22] 
to be quite spectacular, e.g. a deuteron binding energy of Bf ] = 9.6 ± 1.9±i:g MeV in the 
chiral limit, it would appear that their error analysis is incomplete, and when a complete 
error analysis is performed, their claims are likely to be greatly diminished 

V. CONCLUSIONS 

Understanding how nuclei and nuclear interactions depend upon the light-quark masses 
is a fundamental aspect of strong-interaction physics. We have been able to explore the 
m^-dependence of two-nucleon systems using a recently-developed effective field theory and 
naive dimensional analysis. In the ^o-channel we expect that di-nucleon systems, such as the 
di-neutron, are unbound for all values of m q less than their physical values. However, for m q 
larger than their physical values both bound and unbound systems are presently consistent 



WWe should note that EGM apply their NDA argument to the two operators D$ and Dt- However, 
there are four operators, D$\,S2 and T>t\,T2, that can be constructed with one insertion of m q [2]. 
With four operators present, we find that their method saturates the range over which NLO is a 
subleading correction, and thus reproduces the NDA estimates that we have used. 

§§An addendum to Ref. [22] has recently been archived by EGM [24] which increases the error 
in the chiral limit of the deuteron binding energy by a factor ~ 3: B?' = 9.6jZ3 - 2— 2 I— o 4 MeV. 
While their results are becoming more consistent with our analysis, we believe that EGM continue 
to underestimate the errors in the strong interaction parameters. For instance, EGM have not 
implemented the error correlation matrix for parameters in the single- nucleon sector. 
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FIG. 6. The scattering length in the 3 5i-channel as a function of the pion mass for the parameter 
set -1.76 GeV" 2 < d w < -0.91 GeV~ 2 , d 18 = -0.97 GeV~ 2 and rj = 1/20 used in Ref. [22]. 



with data and NDA. In the 3 Si — 3 Di coupled-channels, where the deuteron resides for the 
physical values of the quark masses, the deuteron may or may not be bound in the chiral 
limit. A more definitive statement can only be made with a more precise determination of 
the nN coupling d 16 and a determination of the coefficients of the leading m 9 -dependent 
four-nucleon operators, D 2 . As discussed in Ref. [2], it is likely that a determination of D 2 
will require a future lattice QCD calculation. While we have not investigated higher partial 
waves, given that in this framework they are well described by perturbative pion exchange, 
we do not expect qualitative changes to their behavior away from the physical value of the 
pion mass. It is very exciting indeed to be so close to making fundamental statements about 
nuclear physics. 
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